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Abstract 

We use transport equations to compute the evolution of the quark plasma at the 
electroweak phase transition in the two-doublet Higgs model, obtaining in a more rigor- 
ous and quantitatively accurate way results consistent with previous work. We discuss 
the model in connection with the electroweak baryogenesis scenario, and claim that the 
observed baryon asymmetry of the universe can be obtained within this model with a 
suitable choice of parameters. 



1 Introduction 



Electroweak baryogenesis was proposed years ago |jT[ as a possible mechanism to gener- 
ate the estimated baryon asymmetry of the universe {ub/s ~ 10^^°), and it has been 
the subject of a great amount of work ever since p[. Although the Standard Model at 
the electroweak phase transition satisfies the three necessary conditions pointed out by 
Sakharov 0], a quantitative analysis is far from being trivial. In principle, the actual 
computation of the final baryon asymmetry requires a precise knowledge of the param- 
eters (including those of the Higgs sector), of the details of the phase transition, of the 
i?-violation rate near and at the phase transition, and of the transport properties of the 
QCD plasma in a non-uniform background field, which we do not yet have. A few general 
and persuasive conclusions have nonetheless been reached, and there is some consensus 
on the broad lines of this scenario: 

• At temperatures higher than the EW critical temperature, B violating processes 0] 
are fast, and can wash out any previously created B asymmetry. Soon below the 
critical temperature, B violation is heavily suppressed, and can be safely neglected, 
implying that the baryons we see now must have been created at the EW phase 
transition. This argument assumes the absence of any primordial net B — L number, 
which no EW process (including the anomalous ones) can wipe out and which would 
still be observable today, without any need for baryon production at the weak scale. 
However, in the absence of a compelling model predicting such a. B — L asymmetry, 
we are going to ignore this second scenario. 

• A sufficiently strongly first-order PT is required, to prevent the newly created baryon 
asymmetry from being immediately washed out by the still strong i?-violating pro- 
cesses. 

• During nucleation, the expanding bubbles provide the necessary violation of CP 
and (macroscopic) CPT, thereby creating a net amount of CP-odd fields. This 
process takes place at the bubble walls, where the variation in the background fields 
is most rapid. 
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• The actual B generation takes place when the CP-odd fields reach outside the 
bubble, where 5- violating processes are fast. For this reason, it is important to 
understand the transport properties of those fields. 

• In any case, the Minimal Standard Model is quantitatively inadequate because it 
does not contain enough CP violation so that extensions to the MSM have to 
be considered, and can be constrained (albeit loosely) by this cosmo logical test. 

In this paper, we focus on the transport properties of the plasma at the EW phase 
transition in the two-doublet Higgs model and discuss their effect on baryogenesis. In 
the next section we introduce the model and discuss its qualitative properties near the 
phase transition. In section 3 we show how to derive kinetic equations from microscopic 
principles (a result already present in the literature and apply them to the problem 
at hand. Section 4 contains the explicit computation of the several parameters in the 
kinetic equations, and in section 5 we discuss our results and compare them with previous 
work. 

2 The model 

As mentioned above, we need to extend the Minimal Standard Model to include more 
CP- violating terms. Our choice is to introduce a second Higgs doublet, with the following 
potential |0]: 

V(0l,02) = M4<Pl-viy + X2{4<p2-vif + X3[i<p\<Pl-vl) 

+ i4<p2 - Vj)]' + Mi4<Pl)i4h) - (0l02)(0|0l)] (2.1) 

+ X5[Re{(j)\(j)2) - V1V2 cos^]^ + XQ[Im{(l){(t>2) - V1V2 sin^]^ 

Yukawa interactions couple up-type quarks to 0i, and down-type quarks to either 0i or 02 
(in order to avoid flavour changing neutral currents); it is irrelevant which one we choose, 
since we are going to neglect all couplings but the top's. The extra source of CP violation 
is provided by the angle ^, which cannot be rotated away unless A5 = Xq. 
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What we are really interested in, however, is the effective action, which includes ther- 
mal and quantum corrections. Extremizing the effective action, we can obtain the profile 
of the bubble wall; this computation is discussed in [^] (at least in the case Ai = A2, 
Vi = V2, in which the two neutral Higgses undergo a phase transition simultaneously). 
For the moment, however, we do not need an explicit solution, and write the VEV's in 
the general form: 

< >= pi{x) exp{-te,{x)) (2.2) 

In the following, we shall use h^, {i = 1, 2) to refer to the four (complex) excitations 
in the Higgs fields. 

An important quantity here is the width of the wall compared to the particle mean 
free path; the latter can be estimated using the cross section of the typical QCD process 
at this energy scale, i.e. gluon-mediated quark-quark scattering: as we shall show later, 
Tqcd ~ T/3. On the other hand, the wall width is strongly model-dependent, and even 
in the specific theory we are discussing it has not yet been computed precisely. For this 
reason, it would be useful to develop a formalism able to deal with both the "thin wall" 
and the "thick wall" limits, as well as the intermediate regime. Attempts in this direction 
will be discussed in the last section. Our work here applies, at least in its simplest form, 
to the thick wall limit, where the particles inside the wall can be taken to be in local 
kinetic equilibrium. Under those conditions, it is expedient to follow |^ and perform a 
hypercharge rotation of the fields: 

/. e2^^»^^(^)/, (2.3) 

where /j is a generic field, and yi its hypercharge. 

After the rotation, the spacetime-dependent top mass becomes real (as we said, we 
are neglecting all other masses), but a new interaction term appears: 

Chyp = -2d^ejl (2.4) 

where Jj is the hypercharge current. 

Our task is to study the effect of this term on the plasma density; let's first recapitulate 
the essential features of the system: 
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• The EW phase transition takes place at a temperature Tc ~ lOOGeV^. At this energy 
scale quarks and gluons are deconfined, and the strong coupling constant is small 
enough (a^ ~ 0.1) for perturbation theory to be reliable. 

• The Hubble expansion rate {H ~ T'^/rripi) is much smaller than any other interaction 
rate, and can be ignored. 

• We shall assume that each particle species i is locally in kinetic equilibrium, so that 
its density is determined by its chemical potential fii{x) and temperature. 

• In principle, baryon production depends on the densities of all of the left-handed 
baryons and leptons. However, we shall ignore all leptons and the two lighter quark 
families, since their Yukawa couplings are much smaller than the top's. 

With those assumptions, the system can be described in terms of just three quantum 
number densities: 

q{x) = tL{x) + bL{x) (2.5) 
t{x) = tR{x) (2.6) 
h{x) = hl{x) + h^{x) + hl{x) + h2{x) (2.7) 



Below, we shall derive evolution (i.e., transport) equations for the quantities ( p.5|) - 



(|2.71) , in the presence of the non- uniform external interaction (|2.4| ). Unlike 0, which 
addressed the same problem, we compute the relevant parameters in terms of microscopic 
physics, using a general technique that can be applied to other (possibly more realistic) 
extensions of the Standard Model. 

3 Transport equations 

We refer to ^ for a derivation of macroscopic transport equations from microscopic 
principles; here we limit ourselves to presenting the final result and discussing its physical 
significance. 
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The key idea is that for quasiuniform systems (such as ours), it is possible to define 
local densities r;,j(k, x^) of particles with given on-shell momentum. The time evolution of 
these quantities can be derived by computing appropriate Green functions {i.e. Feynman 
diagrams) involving the elementary fields. 

In the simple case of a scalar field with lagrangian: 

£ = ^d^<j){x)dy{x) - ]^m'<p\x) - V{<P{x)) (3.1) 

we obtain 

(9n(k,x'^) „ „, doj dnCk.x^) , ,,,, 

' + V • Vn k, x'^ + = We{l + n k, x^)) - W^.n k, x^) 3.2 

at ox^ ok'^ 

where u is given by the space-time dependent dispersion relation, and We and Wa are the 

emission and absorption rates: 

Wa{k) = ^Y.\<^\3i^)\^>\'pnn{2rr)H\k-pi+Pr,) (3.3) 
Zuj , 

l,n 

= ^El<^b(0)|/>|'p«K2vr)^5^(A;-pz+p„) (3.4) 

l,n 

Here p is the density matrix of the system, and the sums run over a complete set of states 
of the system. Eq. ( |3l^ ) applies to quasiuniform systems, but does not rely on the coupling 
constants being small, so that in principle we can compute Wa and We to any order in 
perturbation theory, provided that we use the correct finite-temperature Feynman rules 
(for a thorough discussion of the finite-temperature formalism, see [pUj]). 

This procedure can be easily extended to more complicated cases, and to fermions; in 
the latter case, the statistical factors of (1 + n) become {1 — n). For each particle species 
we obtain an equation of the form ( p.2| ) and every scattering process contributes a term 
to the We (Wa) of each emitted (absorbed) field. 

In the case of kinetically thermalized particles, the reaction channel Ai...An Bi...Bm4' 
contributes to the functions Wa and We corresponding to (j) as follows: 
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Zuj J {2Tc}''2hi ^ ^ 



so that the r.h.s. of Eq. (|3.2| ) becomes: 

A,B y'^'^l '^^i B A A B 

(3.8) 

We thus recover the usual detailed balance conditions, and the (local) equihbrium density 
of a given particle species depends only on its dispersion relation. 

Therefore, given a space-time dependent perturbation, the behaviour of the system is 
affected in two ways: 

• Dispersion relations are modified, and that changes the local equilibrium densities; 

• Source terms We and Wa appear, and they determine the rate at which local equi- 
librium is approached. 



In order to find the total number density, we integrate Eq. ( |3.2D over k. Elastic reac- 
tions (which do not change the identity of the involved particles, but only their momenta) 
give rise, after the integration is performed, to the usual diffusive behaviour. The inter- 
actions that do violate some particle quantum number, on the other hand, lead to net 
source terms (they contribute to diffusion as well, but we shall ignore this effect, because 
of its smallness). 

In the next section, we shall apply Eq. ( p. 21) to our problem, considering both effects 
in detail, and presenting some explicit computations. 

4 Results 

The dispersion relation of a particle i is affected by the interaction Chyp as follows: 

PA;^ = ml^ (P + 2yid^ef = m\ (4.1) 

As mentioned above, this automatically provides values for the local equilibrium densi- 
ties i]: 
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f,.„ - ^ (4.2) 

where Ki = (number of spin degrees of freedom) x 2(1) for bosons (fermions), and fli = 
-2%. 

Assuming kinetic equilibrium, the local densities are of the form 

g-/3{£;,{k)-M.) 



where Ei = y + mf — fii. 

In the thick wall approximation, 6 varies slowly, so that we can ignore the third term 
on the l.h.s. of Eq. (|3.2D, which involves the second derivatives of 6, and we only need 



compute the source terms We and Wa- 

• Quark elastic scattering is dominated by QCD processes (see Fig.[^). 



(a) (b) 

Figure 1: Dominant processes responsible for quark diffusion. 

To lowest order, these processes can be computed assuming a perfectly thermalized 
system. In this case we do not have to use Eq. ( p.2| ) and can instead apply standard 
finite-temperature field theory p!0| . 

The finite-temperature rate of quark-scattering processes can be obtained (via the 
optical theorem) from the imaginary part of the quark self-energy. This computation 
is discussed in |]rT|; as emphasized there, diffusion depends not on the total cross 
section / da/dQ, but on the transport cross section | / da/dQsin^ 9, where 9 is the 
scattering angle (qualitatively, this is because close-to-forward scattering does not 
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contribute much to diffusion). Taking that into account, we obtain the interaction 
rate for a particle with momentum p:[| 



where p = |p|, = 4/3, and = (Cy + Nf/2)g^T^/Q, with Cy = 3 and A^/ = 6. 
The leading contribution comes from the infrared behaviour of the integrand, so we 



can write: 

rdiff ( \ 9 2 2t.^ /"^ , ^ 1-cos^g f dk^ 

r,_,(p) - TCjjcose^^j 2(^4 + ^,4 ,,,2^) (4-5) 

and, in the leading logarithm approximation, 

rS{.(p) = Y^C,{Cy + ^)a,2lna;i ~ 0.4- ~ - (4.6) 
(the last step was obtained by averaging over p). 

Higgs elastic scattering is dominated by electroweak processes (Fig.^): 

Diagram 2a is similar to those computed above, except that all of the particles here 
are bosons. A further difference is due to the vector boson mass, which provides 
an IR cutoff: the W and Z bosons receive mass not only from quantum corrections, 
but also from the coupling with the Higgs VEVs. However, both contributions to 
rn^z proportional to Qy/T'^i and in the leading log approximation the precise 
coefficient is irrelevant. The SU{2)l and the U{1)y sectors each contribute with a 
term: 



r%Up) = y'^CACv + ^)aHna-' (4.7) 



i/mS^ in PI is related to We and Wa by: 



r = Weil + ef^^) 
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(a) (b) 

Figure 2: Dominant processes responsible for Higgs diffusion (wiggly lines stand for elec- 
troweak bosons, dashed lines for Higgses, and double lines for any particle coupled to 
electroweak bosons). 

In the SU{2)l term Cy = ^, Cp = 2, a = g'^/in and Nj is the number of doublets 
coupled to the bosons; in the U{1)y term Cy = 0, Cp = I, a = g''^/4:7i and 
Nf = 2J2yi, summed over all the hypercharged particles. The total diffusion width 
from Diagram 2a is 

TniUp) - (4.8) 

The Higgs being a boson, Diagram 2b has apparently the same IR behaviour as 
Diagram 2a, and has to be considered as well; however, in this case the exchanged 
particle has a mass ~ 0(AT), where A is the typical Higgs self-coupling constant 
(see Eq. (|2.1|) ); if A ~ 1, then the diagram contributes a sub leading term, which 
we neglect. 

Both for quarks and for Higgses, the diffusion constant is related to F by the usual 
formula: 



D = (4.9) 
top-quark mass coupling, which mixes q and t (Fig. 3a): 

Strictly speaking, we cannot compute an interaction rate associated to such dia- 
grams, because a quadratic perturbation leads to particle oscillations, not to expo- 
nential decay. However, it is still possible to compute an average rate (which does 
not vanish, since the particles experience thermal damping as well) as follows: 
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As described in [|n[], the top quark retarded propagator in momentum space, 
including thermal correction, has the form: 

i \{pq - So) 7° - p • 7 (l 



P 



S3 
P 



(4.10) 



{po - So) - p2 i^l 

where p = |p| and Sq and S3 are functions of Pq and p. The formula above is 
valid for massless fermions, and does not mix helicities. 

Adding the top mass introduces an off-diagonal term in the helicity basis: 



G 



LL 



G 



RR 



(po-So)7°-p-7(l-f) 



(po-So)'-p2(l-|i) -m 



and 



G 



LR 



G 



—irrit 



RL 



(po-So) -J)M1-^ 



mt 



We now Fourier transform Glr with respect to time and obtain: 



GLR{t,p) 



-Tit/2 



(4.11) 



(4.12) 



(4.13) 



where ul = p'^{l-^^/pf = p^+glT^Cp and = 2/mSo = Us-k'^CfT/Q log(a;i) 
0.5T is the total scattering rate with the thermal bath. 

we can then compute the probability of finding a right top tji at time t + At 
coming from a left top at time t: 



P{R-^L){At) = \GLR{At)\' 



m 



2 



(4.14) 



finally, we can compute the rate of change of the right and left top quantum 
numbers; this is the sum of two terms, since to the time derivative of P{R L) 
we must add the decay products of the left top itself, which have the same 
quantum number. We therefore obtain: 



mt 



-Tht 



(4.15) 
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and, using Eq. (K 



tRip) = (r2 + 4^2) (^^(P) - ^«(P)) (4-16) 



In the limit of large r^, this rate vanishes. As mentioned in |T^, this happens 
because thermal interactions destroy quantum coherence, which is a necessary 
condition for these processes to occur. 

Integrating over p we finally obtain, in terms of variables (7)- (9): 

i = ^^(g - qe, - 2(t - te,)) = r_ (l - ^ - ^) (4.17) 



where Qeq, t^q, and h^q have been defined in Eq. ([4.2[) 



Higgs mass matrix (Fig. 3b): the Higgs mass matrix should in principle be ob- 
tained from the quadratic term in the expansion of the effective action around the 
bubble solution; both outside and deep inside the bubble, the Higgs VEVs are slowly 
varying, and we can obtain the Higgs masses by expanding the finite temperature 
effective potential, Vt{4>) (we assume that the transition is first order, so in both 
regions Vt{4>) is convex and the masses are well defined). In the bubble wall, where 
the field configuration is more strongly space-time dependent, we have to expand 
the effective action. Unfortunately, we do not have an explicit expression for Vt(0), 
nor for the effective action, so we shall just assume that the mass matrix elements 
are of order T. Higgs number violation is due to terms proportional to hihj or hjh'j, 
which are responsible for two point Green functions Ghihy Explicitly, the quadratic 
part of the neutral Q Higgs self-coupling can be written (in Fourier space) as: 

£f^.4^t„)J*^-M ^^^^y^y, ,4,8) 

where 

( h\\ 



H 



hi 



(4.19) 



^The electric charge is not spontaneously broken, so terms of the form do not appear, and we 



can disregard charged Higgses altogether. 
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and M and A are 2 x 2-matrices. 

Thermal interactions modify the dispersion relations by adding "mass-like" terms 
that are much smaller than the original masses, so we neglect this effect; thermal 
damping is taken into account by letting 

kQ^ko = ko + iT/2 (4.20) 

where F is the total Higgs scattering rate computed above. 

A, which is responsible for Higgs number violation, will be taken to be small (in the 
opposite regime, we could content ourselves with the local equilibrium approxima- 
tion); under that assumption, and in the basis that diagonalizes M, the propagators 
are: 

\G,^,f = e-^*sin^ ( ] (4.21) 



2 

mf 



and the same procedure as above yields to: 



" = -rSJL.^ (4.22) 



where 

p(l) _ n^(|All| + |A22|y . 
i Higgs - 4((3)r ^ ' 



Higgs Yukawa couplings (Fig. 3c): 

We shall assume that the Higgs that couples to the top is sufficiently light {rrih < 
mt—rrib) or heavy (m^ > mt+mb), so that either tR — > b^h'^ or h" — > bLtn can occur 
on shell; then, the leading contribution to |jp is simply A|, with a normalization 
factor of 2m j for each fermion. 

In the first case, we obtain: 

2nS (^Eh -Et + \/ml + {kh - ktf^ [{uh + 1)(1 - nb)nt - ^/.^^(l - Ut)] 
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As we said, we assume kinetic equilibrium distributions, so: 



(4.25) 

where and Em are purely determined by the kinematics: 



= E r' - 4W^m V'- ' (4.26) 





— miy 


— Arrifml 


2m? 


Y(m? + m^ 


— mlY 


— Arriiml 


2m? 



EM = E r'^2, '^' + ^m r' - ^ (4.27) 

(in our case, we can safely ignore m^, which is much smaller than the other masses) 
Therefore: 



SXirritmb ( h q t\ , ^ ( h q t\ , 

h = t= -b= * * {--- + -] A = -Tyukawa - - - + - 4.28) 

27r3T \8 6 3; ^ \8 6 3/ ^ ^ 

where 



A{t bh+) = (3 / dE- ^ log ---^—--^ (4.29) 



is of order 1. 

In the second case, Eq. (|428|) still applies, with 

Aih- ^ht)=(3 dE-—-— log (4.30) 
and the roles of irtt and m/i in Em and -Em are interchanged. 

If — > is also allowed (which is true at least in some region of the bubble wall, 
where the Higgs VEV, and consequently m^, are sufficiently small), then we have 
another term, larger than ( ^.30| ) by a factor of nit/ mi,, which therefore dominates. 

In the narrow window m^ — mf, < mh < rrit + mb, all these processes are kinemat- 
ically forbidden and we must go to higher order, including QCD corrections, like 
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the one shown in Fig.|^. The evaluation of this graph proceeds as in the previous 
computations; the exchanged particle is a fermion, so that both its propagator and 
its statistical factor make the diagram less IR divergent; the IR cutoff A is provided 
by whichever is largest among the zero-temperature quark masses, their thermal 
corrections, and |m/j — mj|. In any case, the result is suppressed, with respect to 
the previous case, by a^: 

X^ujpi\mh — mt\ 

^yukawa ~ (4-31) 

Higgs cubic coupling: the effects of cubic Higgs couplings (Fig3d) depend on the 
Higgs masses; if at least one process of the form hi hjh^ can occur on shell, its 
contribution can be evaluated as above (with the small modifications due to the 
different statistics): 



h = V2|A„fep/ , "^^J , , (27r)^,5^(k,+kfc-k,)- (4.32) 

I 7 (27r)32Ei (27r)32Ej- (27r)32Efc^ ; V j t fc ^) v ; 

2Txb (Ej -Ei + ^ml + (ki - kj)2^ [{rij + 1)(1 + nk)ni - njnk{l + rn)] 

here the sum includes all combinations of Higgs and anti-Higgs fields that are kine- 
matically allowed, and Xijk is a function of Ai.-.Ae and of the Higgs VEVs. 

We therefore obtain: 



(4.33) 

The integral, formally divergent as mj 0, is of order T if the Higgs masses are 
~ T, as we expect. 

If no 3-Higgs process is kinematically allowed, electroweak corrections must be in- 
cluded (Fig.|). 
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In this case, 

rgU ^ / J^J [(p + - co*(|) ,4.34) 

Supposing that all the masses are ~ T, so that there is no IR divergence in the 
integral, and remembering that ImU ~ aewT, we obtain that 



2 [-^"feP 

^Higgs ~ Ctew^ijk^^ (4.35) 

with a coefficient of order 1. 
Putting everything together, we obtain the following transport equations: 

Q quark v Q J- yukawa \ q g g I mass \ q g g I (^4.o0J 

I J-^ quark V + J- yukawa \ „ „ o I mass \ „ „ „ I 1^41.0 ( J 

yo 8 3y 8 6 J 

h = Dniggs^'^h + Tyukawa " "g ^ 3^ ~ ^Higgs ^ — ^ — (4.38) 

where Dquark, Dmggs, T yukawa, ^ Higgs = ^^Higgs + ^^Higgs ^re giveu in Eqs. (|]6D, (113), 

(CT) , o), (ii), s). 



5 Discussion and comparison with previous work 

Eqs. ( |4.36| )-( |0^) were ffist derived in 0, who interpreted Chyp as an effective chemical 
potential and used elementary thermodynamics. They made reasonable assumptions to 
estimate the F and D coefficients, as well as both weak and strong sphaleron effects, 
without explicitly computing any of these quantities. In our approach, sphaleron con- 
tributions, which cannot be computed in perturbation theory and have been ignored in 
all our discussion, have still to be introduced by hand. However, they only introduce 
additional terms in Eqs. ( [4.36|) -( ^^38|) , without affecting the terms we did compute. 



One might still wonder whether such a detailed computation is anything more than 
an academic exercise, given that there is no compelling evidence to support this specific 
model, let alone to constrain its many parameters. Nevertheless, in the course of our work 
we have reached a few conclusions that are worth pointing out: 
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Generalizations: the technique we have used can be apphed to any other model; 
admittedly, the thick wall approximation has simplified our task, but it can be 
relaxed: Eq. ( |3.2| ) is valid if the wall is thicker than 1/T, a condition considerably 
weaker than the usual "thick wall" assumption, which requires >> r ~ l/a^T. If 
the latter constraint is not satisfied, but the first one is, we can still use Eq. ( fj.2D , but 
we cannot assume local kinetic equilibrium. This means that the local densities are 
not simply parametrized by chemical potentials and the whole computation 

becomes quite cumbersome, but hopefully still tractable. 



We should mention here that ||T2| tried to avoid this complication, deriving a unified 
formalism that does not require the wall to be thin or thick. They derived transport 
equations similar to ours, by adding source terms (corresponding to our heq,teq,qeq 
terms) to the unperturbed evolution. The sources were computed for generic wall 
thickness, but taking only quadratic interactions into account; for this reason, they 
found vanishing effects in the thick wall limit, at odds with our own result. For 
realistic values of Ly^ the numerical disagreement may not be important, but still it 
is not negligible. 

Diffusion: The quark and Higgs diffusion constants that we obtained are noticeably 



smaller than the estimates in |13]: Dq^ark = 6/T, Dmggs = 110/T, whereas they are 
comparable with the values used by for their numerical analysis: D quark = 3/T, 
D Higgs = 10/T; for this reason we believe that the qualitative results and the 
discussion contained in are still valid. 

Source terms: in 0, Twiggs and Tmass were estimated to be: 

^ Higgs ~ ^mass ~ ^ (5.1) 

As we said earlier, Tmass is sensitive to the quark mean free path, and the value 
we obtained in Eq. ( |4.17| ) is considerably smaller. Twiggs strongly depends on the 
unknown parameters of the Higgs sector; Eq. ( [4.33| ) can be used to study the whole 
parameter space, but realistically we won't obtain results that are too different from 



16 



• Baryon number: ^ have integrated Eq. ( |4.36| )-( P.38| ) numerically, and found that 
values of ub/s consistent with experimental estimates can be obtained in this model 
by taking ^ ~ 10^^ — 10^^. We have not performed the same calculation with our 
values of the F's and the D's, but we believe that their conclusions would not be 
qualitatively altered. Therefore, this model should still be regarded as a viable 
candidate for electroweak baryogenesis. 
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Figure 3: Vertices that violate at least one among quantum numbers t, g, h. 
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